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Investigations of Bragg reflectors in nanowire lasers 
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The reflectivity of various Bragg reflectors in connection to waveguide structures, including nanowires, has 
been investigated using modal reflection and transmission matrices. A semi-analytical model was applied 
yielding increased understanding of the diffraction effects present in such gratings. Planar waveguides and 
nanowire lasers are considered in particular. Two geometries are compared; Bragg reflectors within the 
waveguides are shown to have significant advantages compared to Bragg reflectors in the substrate, when 
diffraction effects are significant. 
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I. INTRODUCTION 

Semiconductor nanowires including nanowire lasers are 
promising as building blocks for realization of nanoscale 
photonic devices^. Various techniques such as molec- 
ular beam epitaxy (MBE) or metalorganic chemical 
vapour deposition (MOCVD) can be used to form 
nanowires with accurately controlled geometry and ma- 
terial composition, yielding a high level of flexibility— 
To obtain an efficient laser resonator, the reflectivity at 
the end facet of the nanowire must be high. The refrac- 
tive index contrast between the semiconductor nanowire 
and its surroundings is typically very large; the simplest 
nanowire-laser designs could thus use the cleaved end 
facets as reflectors. With such designs, the reflectivity 
of the guided lasing mode is quite moderate for single 
mode semiconductor nanowire lasers (~ 25% for GaAs- 
based nanowire, ~ 18% for ZnO based)^. Bragg grat- 
ings have been proposed to obtain a higher end facet re- 
flectivity. Such gratings are fully compatible with most 
nanowire fabrication methods, e.g. MBE and MOCVD, 
and have already been realized experimentally^. Chen et 
al. have performed numerical analyses of nanowire Bragg 
structures. They show that a nanowire supcrlattice can 
be used to achieve near unity modal reflectivity at sin- 
gle mode operation^. Additionally, they have performed 
an optoelectronic analysis of nanowire lasers with dis- 
tributed Bragg reflector mirrors^, showing a significant 
improvement in output power. Fricdler et alr^ have used 
coupled mode analysis to calculate the reflectivity of a 
dielectric Bragg grating within a GaAs nanowire. They 
conclude that the reflectivity of such a grating is rather 
poor in the single mode regime for a GaAs nanowire, and 
propose to rather use metallic mirrors. 

When the lateral scale of a waveguide is of the order 
of the wavelength of the guided light, diffraction effects 
become significant. In this work we perform a detailed 
analysis of Bragg grating reflectors in connection to such 
diffractive waveguides, to investigate in which regimes a 
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FIG. 1. Two geometries for implementation of Bragg mirror 
at the end facets of a nanowire laser. 



Bragg grating is efficient. A semi- analytical model will 
be used; as compared to finite clement methods this helps 
in explaining more of the mechanisms that influence the 
reflectivity. A substrate grating (Fig. [T] II) is found to 
have surprisingly low reflectivity compared to a grating 
within the waveguide (Fig. [T]I). Furthermore we see that 
even for extremely small waveguides, where only a small 
fraction of the field is within the waveguide, a near unity 
reflectivity may be obtained by having enough periods in 
the Bragg grating. 

The reflection and transmission properties of an inter- 
face are fully described by reflection and transmission 
matrices. These matrices describe the amount of mode 
i that is reflected or transmitted into mode j. In our 
previous work a formalism was developed to calculate 
the transmission and reflection matrices for end facets 
of waveguides^. The method is particularly useful for 
highly diffractive waveguides and nanowire laser applica- 
tions. In this paper we extend the formalism to examine 
the effect of Bragg reflectors, and consider two differ- 
ent geometries; (I) Bragg reflector in the substrate, (II) 
Bragg reflector at the top of the nanowire (within the 
waveguide). The two geometries are sketched in Fig. [T] 
A heterostructure based on GaAs and Alo.3Gao.7As is 
used throughout as an example. 

The outline of this article is as follows: A presentation 
of the multimodc transfer matrix formalism is presented 
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in Sec. mi The calculation model for the reflection and 
transmission matrices is presented in Sec. IIIII A brief 
summary of previous work describing the reflection at 
the end facet of a waveguide is given, as well as the nec- 
essary generalizations to describe a Bragg grating within 
a waveguide or in the substrate. Sec. IIVI contains a dis- 
cussion concerning the design of the Bragg gratings. Nu- 
merical results for a planar waveguide structure arc given 
in Sec. El and some results concerning nanowircs, with 
2D confinement are given in Sec. IVII 
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II. TRANSFER MATRIX FORMALISM 

The theory of transfer and scattering matrices can be 
found in standard textbooks^. We will here briefly re- 
view the concepts, to introduce our choice of notation. 
Consider a stack of layers, with layer boundaries perpen- 
dicular to the propagation axis, z. Each layer is homoge- 
neous w.r.t. z. The field in each layer can be described 
using its modes. Throughout this article we define modes 
as being pairs of electric and magnetic fields that are 
cigenfunctions of the electromagnetic propagation oper- 
ator along the z-axis. For an infinitely long waveguide 
that is homogeneous along the z-axis, the modes will 
correspond to the eigenmodes of the whole structure. 
However, for waveguides of finite length or with inho- 
mogeneities the modes are merely local modes, not to be 
confused with the supermodes of the overall structure. 

Let the forward propagating mode n in layer i have am- 
plitude a % n , and the backward propagating mode have am- 
plitude b l n . The vectors a 1 and h 1 contain the amplitudes 
of all forward propagating modes and backward propa- 
gating modes, respectively. Let r 31 and t 31 be matrices 
describing the modal reflection and transmission respec- 
tively, for light incident from layer i towards layer j; simi- 
larly r 13 and t 1J describe the reflection/transmission from 
the opposite side. Using these matrices, we can relate the 
field in layer i to the field in layer j: 



r%4 + J2l%4 

v 
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FIG. 2. The transfer matrix M Jl relates the field in layer i to 
the field in layer j. 



plicitly expressed using the field in layer i, i.e., 



a 1 



(4) 



The matrix M 31 is known as the transfer matrix; a gen- 
eral transfer matrix is illustrated in Fig. [5] In light of (Q]) 
it can be expressed in terms of the scattering matrix S J% : 



M 31 



/•" /•'< (/'••') ' r". r ij (/'••') ' 

-(t ij ) + r ji , (t ij y 



(5) 



Here, the superscript (+) denotes the matrix inverse 
or More-Penrose pseudo inverse depending on whether 
the matrix is quadratic or rectangular. The presence of 
evanescent modes in a layer could cause ill-conditioned 
transfer matrices. As the transmission coefficients t 1 ^ in- 
volving evanescent modes may be extremely small, ma- 
trix inversion of the transmission matrix may cause nu- 
merical instabilities. To avoid such problems it is prefer- 
able to use recursive relations derived from transfer ma- 
trices rather than direct matrix multiplication. We con- 
sider a stack of three layers, 1,2 and 3; the individual 
transfer matrices are multiplied to find the total reflec- 
tion and transmission properties. Recall that r Jl ( t Jl ) 
denote the reflection (transmission) from layer i to layer 
j. The combined reflection and transmission coefficients 
for the system of layers are given by: 



rC, i.e. r 3% = r"y , similarly 

for t 3% , r lJ and t %3 . Eq. |T]) can be rewritten in matrix form 
as 



= S 31 



Here S 31 is the scattering matrix: 



r 3 ' 1 t 13 
pi r ij 



(2) 



(3) 



When considering a sequence of layers, it is convenient 
to reformulate © so that the field in layer j can be ex- 
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r 3 V 2 ) r 3 ^ 1 (6a) 

r-r 32 r 12 )V 2 )* 21 (6b) 

I - r 32 r 12 ) + t 23 (6c) 

l ) + t 23 . (6d) 



Propagation in the z-direction within one layer can 
be described in the same manner. Mode k propagates 
according to 



a" u = e ^" a 



k ' 



(7a) 
(7b) 
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FIG. 3. Circular waveguide of diameter 2p oriented along the 
z-axis with an end facet at z = 0. The waveguide core and 
cladding have refractive indices n co and n c i, respectively. For 
z>0 there is an ambient medium, which is assumed to be ho- 
mogeneous with refractive index ra a - The aim is to determine 
the reflection and transmission matrices for the end facet, as 
illustrated in the figure. 



where fik is the modal propagation constant in z- 
direction, and d is the propagation distance. 



III. FINDING THE REFLECTION AND TRANSMISSION 
MATRICES 

The problem of finding the reflection and transmis- 
sion matrices at the end facet of a waveguide terminated 
in a homogeneous medium has been addressed by us 
previously, here we briefly sum up the main results. The 
geometry of the problem is shown in Fig. [31 here exempli- 
fied using a circular waveguide. We describe the field at 
both sides of the boundary, z = 0, using a set of modes. 
The modes in the half-space z>0, constitute a continu- 
ous set of radiation modes, whereas for z<0 the modal 
spectrum consists of a discrete set of bound modes and a 
continuous set of radiation modes. The modal spectrum 
is discretized using periodic boundary conditions at each 
side, the width of the computational cell in both x and 
y direction is 2L. The electric field of mode m in the 
ambient half space can be written as 



£ m {x, y) = £ m exp{ik x x + ik y y). 



(8) 



The magnetic field, % m , is described in the same way. 
The label m is a collection of the modal indices, m = 
(p, g,pol). The polarization, pol, is TE or TM, and the 
real transverse wavevectors are k x 



L ' 



where 



is given by fc 2 = rfyj 1 /c 2 — fc 2 — fc 2 , where n a is the refrac- 
tive index of the half-space z > 0, and c is the vacuum 
light velocity. The constant vectors can be expressed^ 2 - 



£m(TE) = Ay/UfI{-ky, k x ,0), 

A 



H 



i(tm) = yk x k z ,k y k Z) - (kl + fc 2 )) 

m(TE) = —=(-k x k z , ~k y k z , kl + fc 2 ), 



^m(TM) = A^/uJs2(-ky, k x ,0), 



(9a) 
(9b) 

(9c) 
(9d) 



where A = 1 j y(fc 2 + fc 2 )|fc 3 |2L 2 . We assume that 

the medium is nonmagnetic, and the permittivity of the 
medium is e a . The modal fields of the waveguide are de- 
noted e< = ei(x,y) and h; = hi(x,y), i = 1,2, . . .. We 
now use the continuity of the transverse electric and mag- 
netic fields. Assuming the incoming mode {ej,hj}, we 
can write 



a,wg (t) 



E,a,wgc(t) 
mi 



(10a) 



hf -E^ WSh ? = Ec g ^> ( 10b ) 



valid for all x and y. Here r^ wg is the reflection coeffi- 
cient from mode i to mode j, and t^ s is the transmis- 
sion coefficient from mode i (z < 0) to mode m (z > 0). 
The superscript (t) stands for the transverse component 
(x and y components) of the vector. Eqs. (|10a[) and 
(|10b[) can be combined as follows. Take the vector prod- 
uct between (|10a.p and H^*(x, y), and integrate over the 
unit cell. Similarly, take the vector product between 

£^,*(x,y) and (|10b[) . and integrate over the unit cell. 
Combining the resulting equations yield 

r a ' ws = b,] a ' Wg = [$im - tt^HSfa, + * im ] + , (Ha) 

t a ' Wg = [t mi r WS = \ ([$mi + - [*m< - *mi]r a ' Wg ) . 

(lib) 

Here we have defined the inner products 
K(m)* 1 



|«(m)| 2 



J ef xU%*{x,y)-Z&A, (12a) 



cell 



= / *) x hf* • zdA (12b) 

|/c(m)| 2 7 

cell 

The unit vector in the z-direction is z, and 

< m ) = { kl'pol^TM. ( 13 ) 



p and q are integers. The modal propagation constant k z 



It is straightforward to extend the formalism to de- 
scribe the reflection and transmission for light incident 
onto the facet from the ambient medium. We assume 
the incoming wave {£i,Hi}, and consider the boundary 
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conditions, similarly to pop . The resulting expressions 
for the reflection and transmission arc 



t WS ' a = [* mi ] Wg ' a = ^ ([*mi + *,».]' - [*mi - $™]*r WS ' a ) 



(14a) 
(14b) 



Here, 



|/3 m |2 



e^xf (^).Jdi, (15a) 



cell 



/?m 1 / p (t) 



|/3m| 2 



J £®(x,y)xhM*.idA. (15b) 



We have here assumed that the waveguide modes are 
orthogonal, and fulfill 



If B J * 
eP x h J * • zeL4 = f— rrJ 

I \P 3 
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(16) 



t a ' wg = 2 f [$ mi ]+ + [* mi 



This ort honor mality relation can always be fulfilled for 
the modes of nonabsorbing waveguides^. Note that the 
transmission matrices may also be found directly from 
the inner products; 



(17a) 
(17b) 



A mode is said to be real if it has a real-valued propa- 
gation constant and the transverse electric and magnetic 
field of the mode can be written real for all values of x and 
y. For modes in nonabsorbing waveguides the transverse 
fields can always be written reali^. Modes with a real- 
valued propagation constant will therefore be real modes. 
For coupling between real modes i and m with real 
propagation constants, we have = K(m)*/|/t(m)|. 

In this case $ m! and are both real, and we have 
^mi = ^>mi, and $ m i = <t m i. We then see directly that 
^a,wg _ £wg,a^ cxac tly as predicted by the reciprocity the- 
orem. 

It is a necessary condition when solving for the reflec- 
tion matrices to have a well defined system of equations. 
A minimum requirement is to use the same number of 
orthogonal modes at both sides of the boundary. This is 
however not an ideal solution, as the sampling in the spa- 
tial frequency domain is quite different at the two sides 
of the boundary; thus a large number of modes would 
be necessary for an accurate description of the interface. 
We have rather chosen to use a higher number of modes 
on the ambient side; this enables a good description of 



the forward reflection 



vg . However, for the backward 



reflection r wg,a , we cannot directly find the reflection co- 
efficients for all modes of the ambient medium. The pro- 
cedure is as follows: First we find the reflection matrix 



for the ambient modes with the lowest spatial frequen- 
cies. These modes must be sufficiently well described 
by a superposition of waveguide modes. More precisely 

they obey (k% + k%) < ((n co w/c) 2 - f3? im \ , where ftim 

is the possibly imaginary propagation constant of the 
highest order waveguide mode. For higher spatial fre- 
quencies, we approximate the reflection coefficients using 
the scalar Fresnel equations for reflection at an interface 
with index contrast n a /n c i. This approximation shows 
very good agreement provided the number of modes on 
the waveguide side is not too small. 

For Bragg gratings within the waveguide one also needs 
to describe the reflection and transmission properties 
when there is a waveguide at both sides of the inter- 
face. One possibility for performing this calculation 
would be to repeat the procedure described previously, 
using waveguide modes at each side of the interface. As 
we have waveguide modes at both sides of the bound- 
ary the inner products similar to (|12|) and p5p could no 
longer be formulated as Fourier transforms. To avoid 
this problem, we formulate the reflection and transmis- 
sion matrices for transitions between waveguides in terms 
of the previously acquired relations for a transition from 
a waveguide to an ambient pip . We start by formu- 
lating the boundary condition as in (jXDJ . The fields at 
both sides are then expressed in terms of the inner prod- 
ucts (|12p between each waveguide and a dummy ambi- 
ent layer. Using this procedure we obtain expressions 
for the reflection and transmission between waveguides 
formulated in terms of their reflection and transmission 
matrices towards a dummy ambient medium. Note that 
there are no assumptions made here, and the accuracy is 
given from the accuracy of the reflection and transmis- 
sion matrices from the waveguide to the dummy ambi- 
ent medium. The details concerning this calculation are 
given in Appendix [X] For a transition from waveguide b 
towards waveguide c the result is 



t cb = Qcb (J _ r ab r cb^ ( lga ) 
r cb = (peb _ r ac G cb r ab^+ (Qcb^b _ ^CQob^ > ( lgb ) 

where the superscript a denotes the dummy ambient me- 
dia, and 



G cb = (J _ r ac r acj ^ac) + f ab ^ _ ^6^6) + 



(19a) 



The opposite transition is described by interchanging 
indices b and c. 

For a thin diffractive waveguide, the imposed bound- 
ary conditions cause artificial reflections from the bound- 
ary. To deal with these artificial reflections, we introduce 
some loss into the system, i.e., e — > e + ijeq, at both sides 
of the boundary. The loss parameter 7 should be small 
enough not to alter the reflection properties of the bound- 
ary significantly^. Note that this loss is merely artificial, 
and it will only be included when necessary. To describe 
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the Bragg grating, we must treat the scattering at inter- 
faces as well as propagation in homogeneous layers. The 
loss is not included in the propagation description, as 
this would lead to an underestimate of the reflection as 
compared to the physical situation. Since the loss is in- 
cluded in the description of the interfaces, but not in the 
propagation description, it represents a deviation from 
a physical structure and some error is to be expected in 
the final result. Decreasing 7 will decrease this error. We 
previously assumed that the waveguide modes fulfilled 
the orthonormality relation Q16p . This is however only 
generally true for nonabsorbing waveguides. For slightly 
absorbing waveguides we may assume that the deviation 
from (TTB|) is smali^. The orhonormality relation can even 
be exactly fulfilled for the planar step index waveguides 
considered in this paper when e — > e + 1760 ■ 

IV. DESIGNING THE BRAGG GRATING 

In this section we consider the design of the Bragg grat- 
ing structure in connection to a waveguide. Let a waveg- 
uide be terminated by some grating consisting of layers 
with refractive indices n\, and n c , and thicknesses db and 
d c respectively. The waveguide itself has refractive in- 
dex n c . The structure is designed to be a quarter wave 
stack for the fundamental mode of the waveguide, i.e. the 
mode with propagation constant /3c^ in z-direction. The 
thicknesses of the quarter wave layers are given by 

dc = -^T ( 20a ) 
d h = -^ (20b) 

The response of the grating is highly dependent on 
which material that constitutes the terminating layer. 
To illustrate this we consider a planar structure with 
the grating within the waveguide (Fig. [1] I). The het- 
erostructure is based on GaAs and Alo.3Gao.7As and the 
ambient medium is vacuum. The waveguide where the 
lasing is to occur consists of GaAs, i.e. n c = n(GaAs), 
rib = n (AlGaAs). The lasing wavelength for GaAs in 
the Zinc blende (ZB) crystal phase is 870 nm at room 
temperatureM. At this wavelength, the refractive in- 
dices of GaAs and Alo.3Gao.7As arei£ n (GaAs) =3.6 and 
n (AlGaAs) =3.4. The structure is surrounded by air. 
The total reflection matrix is found using the recursive 
relations (JB]) • To this end we need the propagation matri- 
ces, the reflection and transmission matrices describing 
the interfaces between the waveguide layers, and the re- 
flection matrix for the transition from the terminating 
grating layer towards the surrounding ambient. 

We calculate the total reflection when the terminating 
layer consists of either the low index or the high index 
material. Let the frequency of the light be u> and the 
width of the waveguide be 2a. We include all modes 
with P 2 > j3f im , where /3u m is the cut-off limit. In this 



example, a(ui/c) = 1, /3i im = 10i(ui/c), L — 100a, and 
7i =0.1. Fig. 0] shows the reflection coefficient for the 
fundamental TE even mode, as a function of the num- 
ber of periods, for gratings terminated by either GaAs or 
Alo.3Gao.7As. Note that the behavior is fundamentally 
different depending on whether the grating is terminated 
by the high index material or by the low index material. 
When the material with the lowest refractive index ter- 
minates the grating, the reflection is reduced rather than 
increased for the first layers. 



This can be explained as follows. For a quarter wave 
stack all multiple reflections interfere constructively, as 
there is an additional phase shift of 7r at every second in- 
terface when the refractive index goes from low to high. 
If the interface towards the ambient layer breaks this pe- 
riodicity the portion of the field reflected at this last 
interface will interfere destructively with the rest of the 
field. Unlike in conventional Bragg gratings, this last re- 
flection may be crucial, as there is such a large index 
contrast between the grating and the surrounding air. If 
the grating consists of several periods, most of the field 
is reflected before it reaches the last interface; the effect 
of this additional phase shift is therefore gradually re- 
duced. Fig. [5] clearly shows that to enable an efficient 
grating the terminating layer should be made from the 
high index material. An alternative solution if one needs 
to terminate the grating using the lowest index material, 
is to grow the terminating layer with twice the thickness 
to compensate for the phase shift. A waveguide grating 
terminated by the highest-index material followed by an 
ambient of the same material would yield a similar ef- 
fect, as the effective refractive index of the fundamental 
mode in the waveguide is lower than the refractive index 
of the bulk material. This may be part of the reason why 
Fricdler et al^ obtain so low reflectivity for the thinnest 
waveguides. As the thickness of the waveguide increases, 
the index contrast and thus the reflectivity at the last 
interface will decrease, and this effect would diminish. 



The phase shift at the first interface has a similar ef- 
fect. If the first interface breaks the periodicity, the con- 
tribution from this first reflection will be out of phase 
with the remaining contributions. This situation may e.g. 
occur when a GaAs waveguide is terminated by an Al- 
GaAs/GaAs Bragg grating in the substrate (Fig.[T]II). As 
the waveguide thickness decreases below a certain limit, 
the effective refractive index of the fundamental mode in 
the waveguide will decrease below that of the first layer 
of the substrate (consisting of AlGaAs). This will cause 
the reflection at the first interface to interfere destruc- 
tively with the remaining backscattered field, leading to 
reduced reflectivity One possible solution to compensate 
for this phase shift, is to adjust the thickness of the first 
layer accordingly. 
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FIG. 4. Reflection coefficient |r^ tal | 2 of a GaAs/Alo. 3 Gao. 7 As 
Bragg grating in a planar waveguide, with heterostructure 
consisting of media with n = 3.4 and n = 3.6. The structure 
is surrounded by vacuum. The figure shows the reflection co- 
efficient when either the low index material or the high index 
material is the final layer of the structure. Here, a{uj/c) — 1. 



V. PLANAR WAVEGUIDE STRUCTURE 

We now look more closely into some numerical exam- 
ples for a planar waveguide with a Bragg grating. Two 
situations are considered; the grating is either within the 
waveguide (Fig.[TJI) or in the substrate below the waveg- 
uide (Fig.[TJII)). We also briefly consider an intermediate 
geometry. The planar waveguide with ID confinement 
is less computationally demanding compared to the 2D 
case; in addition both bound and unbound modes can 
be found analytically^. The planar case is therefore well 
suited to test qualitative relations and convergence crite- 
ria. In a planar waveguide there is no coupling between 
modes of different parity (odd/even) or between modes 
of orthogonal polarization (TE/TM). The discussion is 
therefore limited to even TE-polarized modes. 

First, let the Bragg grating be within the waveguide, 
as shown in Fig.[TJI). Such structures can be realized by 
growing the Bragg grating at the end of the nanowirc 
growth, by alternating the source materials during the 
epitaxial growth. In this example, the main part of the 
waveguide consists of GaAs, and alternating layers of Al- 
GaAs and GaAs are grown at the top of the waveguide. 
The uppermost layer consists of GaAs, and the structure 
is surrounded by vacuum. We have performed calcula- 
tions for up to 100 periods, to see the behavior in the 
limit of several periods. Note however that this is a very 
high number, which is not easily achieved with today's 
technology. 

We consider four normalized waveguide widths; 
a(uj/c) = 0.1, a(u)/c) = 0.5, a(ui/c) = 1, and a(ui/c) = 
10. For reference, the single mode regime for even 
TE modes in this GaAs waveguide extends up to a = 



Grating within waveguide 
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FIG. 5. The reflection coefficient |r 1 i tal | 2 as a function of the 
number of periods in the GaAs/Alo.3Gao.7As Bragg grating, 
for four waveguides of various width, a. 



0.91(w/c). The computational cell half-width L is 
100/(w/c), and the loss parameter is 7 = 0.1. The 
cut-off limits Aim were chosen by considering the con- 
vergence of the reflection in each case. In the order of 
increasing waveguide width we used Aim = 20(ui/c)i, 
Aim = 15(o;/c)i, Aim = 10(w/c)i, and Aim = 5(w/c)z. 
The reflection coefficient |r'° tal | , i.e. the amount of the 
fundamental mode reflected back into itself, is shown in 
Fig.[5]as a function of the number of periods in the Bragg 
grating. 

It is seen that the end facet reflectivity of all waveg- 
uides converge towards a value very close to 1. Even as 
the normalized waveguide width decreases below auj/c = 
0.5, one can still obtain high reflectivity by increasing the 
number of periods in the grating. As will be seen later 
this is contrary to what is observed for the case with the 
grating in the substrate. 

Before we proceed it is instructive to review some- 
what how the modal fields are influenced by diffraction. 
Firstly, as the width of the waveguide decreases, a de- 
creasing proportion of the modal field will be confined 
within the core of the waveguide. As a consequence, the 
effective refractive index of the fundamental mode de- 
creases towards the limit where it is close to the refractive 
index of the cladding material. Secondly, as the modes 
are confined to smaller areas in space, a corresponding 
spreading of the spatial frequencies of the mode must 
follow. This will e.g. imply that the waveguide modes 
will couple more strongly to each other upon reflection, 
and that there will be a larger angular spread of the beam 
upon transmission towards an ambient medium. 

A grating within the waveguide with a relatively low 
refractive index contrast will roughly preserve the same 
set of modes along the grating. Except for the modes that 
are very close to their cut-off, each mode will thus expe- 



7 



Grating within waveguide 









r 












C) = 1 








a(ct>/ 

































5 10 15 20 25 30 

Number of periods 

FIG. 6. The reflectivity |rl° tal | 2 as a function of the number 
of periods in the Bragg grating for a waveguide of normal- 
ized width a/(u)/c) = 1. The Bragg grating consists of two 
materials with refractive indices n = 4 and n = 2. 

rience a jump in the effective refractive index in a man- 
ner quite similar to what is seen for conventional Bragg 
gratings. Only a small amount of energy will therefore 
be transferred from e.g. the fundamental mode to the 
higher order modes. If the index contrast of the grat- 
ing is larger, each transition in the grating represents 
a more significant perturbation to the modal field, and 
there will be a larger amount of cross coupling between 
modes. Some of the energy from the fundamental mode 
may thus couple into other modes. Fig. [5] displays the 
reflectivity of a Bragg grating consisting of two mate- 
rials with higher index contrast. Here GaAs has been 
replaced by a material with refractive index 4, and Al- 
GaAs has been replaced by a medium of refractive index 
2. The simulation parameters are the same as for the 
corresponding GaAs/Alo.3Gao.7As structure. Note that 
the reflection coefficient of the fundamental mode now 
converge towards a value less than unity. As can be seen 
by comparing Fig. [5] and Fig. ® there is a trade-off here 
in terms of the index contrast. Higher index contrast en- 
ables quite high reflection using fewer periods. On the 
other hand the maximum obtainable reflection is larger 
for the lower index contrast system. 

We proceed to consider a Bragg grating in the sub- 
strate below the waveguide, as shown in Fig. |T] II). For 
nanowirc applications, such structures can be realized by 
growing the substrate Bragg grating before the nanowire. 
This geometry has the advantage that it is easier to con- 
trol the thickness and composition of the layers compared 
to the structure with the grating within the nanowire. 
The structure consists of the same materials as for the 
case with the grating within the waveguide, and we con- 
sider the same four waveguide widths as before. The 
half-width of the computational cell, L, used in the cal- 
culations was 1000/ (cj/c), which is larger than for the 
case with the grating within the waveguide. The rea- 



son for this is that the transition towards this substrate 
Bragg grating represents a more significant change in the 
modal fields. The coupling from the fundamental mode 
to higher order modes including radiation modes is there- 
fore enhanced, and these modes will be more influenced 
by the artificial boundary conditions. The modal cut-off 
limit was taken to be f3n m = 3(w/c)i. The resulting re- 
flectivity as a function of the number of periods is shown 
in Fig. El 

In the geometric optics limit, i.e. as the normalized 
width of the waveguide increases, this substrate grat- 
ing and a corresponding grating within the waveguide 
should approach each other. In this limit the reflection 
and transmission coefficients of the bound modes can be 
approximated by those of plane waves at a homogeneous 
interface^. Comparing Fig. [7] and Fig. we see that 
this approximation is indeed accurate for a = 10(w/c). 
However in the highly diffractive regime, there are large 
differences between the two Bragg geometries. Using a 
substrate grating, increasing diffraction will lead to de- 
creased reflection. It is not possible to compensate for 
this by adding more layers. Fig. [5] helps us understand 
this effect. The upper plot displays the reflectivity of 
the quarter wave stack separately, i.e. the reflectivity 
of plane waves incident from the first substrate layer to- 
wards the remaining quarter wave stack. The lower plot 
displays the transmission coefficients from the fundamen- 
tal mode of the waveguide into each of these plane waves. 
The Bragg grating has reduced reflectivity in the region 
k x = 0.75(w/c) to k x = {uj/c). Increasing the number of 
periods in the Bragg grating will increase the frequency 
of the oscillations in this region, but it will not decrease 
the width of this region with reduced reflectivity. The 
energy transmitted into plane waves in this low reflectiv- 
ity region will therefore be partly transmitted through 
the grating and transported away. This explains why we 
do not achieve high reflectivity. For highly diffractive 
waveguides, a significant amount of the energy is trans- 
mitted into evanescent modes (k x > oj/c). The evanes- 
cent modes do not transport energy, so eventually they 
will couple back into propagating plane waves; especially 
to the plane waves with similar spatial frequencies. It is 
thus natural to assume that most of the energy in the 
evanescent modes is coupled back into the plane waves 
with reduced reflectivity; thus a large part is transported 
away from the structure. As au>/c — > 0, the portion of 
the fundamental mode transmitted into the region with 
k x < 0. 75 (u;/c) will decrease, thus reducing the effect of 
the grating. 

For very small aw/c, the effective refractive index of 
the fundamental mode tends to the refractive index of 
the cladding medium (vacuum). The reflectivity of the 
fundamental mode will therefore converge to zero for 
very thin waveguides if the cladding material is the same 
as the substrate, as would be the case in the absence 
of a substrate grating. By adding a substrate layer of 
Alo.3Gao.7As at the end of a GaAs waveguide, the effec- 
tive index contrast for the fundamental mode will first 
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FIG. 7. The reflectivity |rfi tal | as a function of the number 
of periods in the Alo.3Gao.7As/GaAs Bragg grating, for four 
waveguides of various width, a. 
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FIG. 8. Transmission coefficients from the fundamental mode 
into the plane wave components in the first layer of the sub- 
strate grating (lower plot), and the respective reflectivity 
when these plane waves are reflected due to the grating (upper 
plot). Note that the reflectivity is higher than unity for spa- 
tial frequencies corresponding to propagating modes in GaAs, 
but evanescent modes in the Alo.3Gao.7As. As no energy is 
transported by the evanescent modes, this does not violate 
energy conservation. 



decrease and then increase again as the waveguide width 
is decreased. As a consequence, the reflection coefficient 
for the fundamental mode towards the Alo.3Gao.7As sub- 
strate will vary correspondingly. One may thus achieve 
relatively high reflectivity, but this is due to the fact that 
there is a high effective index contrast at the waveg- 
uide/substrate interface, not due to constructive inter- 
ference in the quarter wave grating. 

As a function of the number of periods, the reflectiv- 
ity for the thinner waveguides decreases before it starts 
to increase (Fig. [7]). This can be understood in terms 
of the phase shifts at the interfaces. As the effective re- 
fractive index of the waveguide decrease, the phase of 
the reflection at the interface between the waveguide and 
the substrate will change. As discussed in Sec. IIVI this 
phase shift may lead to destructive interference between 
the backscattered contributions. For the two extreme 
waveguide widths a = 10(a;/c) and a = 0.1(w/c), the 
phase shift of the first reflection coefficient is close to 
or 7r, respectively This leads to destructive interference 
for the waveguide width a = 0.1(w/c). To compensate 
for this, we doubled the thickness of the first layer in the 
Bragg grating, which strongly increased the reflectivity. 
The result is shown in Fig. [7] 

We have seen that the differences between the two 
grating geometries become large in the diffraction limit. 
Before we proceed to 2D calculations on nanowires, we 
therefore consider an intermediate geometry. Here, the 
segments of the grating have a larger lateral width than 
the central waveguide. Such structures can be realized 



by first growing a substrate Bragg grating, and then etch 
to reduce the lateral size of the reflector. This could be a 
potential way to overcome some of the weaknesses asso- 
ciated with substrate Bragg gratings, while maintaining 
a structure that is relatively easy to fabricate. The width 
of the central waveguide in the calculation was taken to 
be a{uj/c) = 1. Fig. displays the reflection coefficient 
| r totai J ag a f unc tj on f t ne number of periods for a grat- 
ing of lateral width a(w/c) = 2. The simulation param- 
eters were /3h m = 5(w/c)i, ji = 0.1, and L = 100/(o;/c). 
As a reference, we also show the corresponding reflec- 
tion coefficient for a substrate grating and a grating with 
equal width as the central waveguide (taken from Fig. [5] 
and Fig. [7J . 

From Fig.O we see that the reflection of the a{uj/c) = 2 
grating is intermediate between the two geometries dis- 
cussed earlier. There are also large fluctuations as a func- 
tion of the number of periods. In the transition from 
the central waveguide (a(w/c) = 1) towards the segment 
with a(uj/c) = 2, the field of the fundamental mode ex- 
periences a significant alteration. This leads to large cou- 
pling into several higher order modes of the a{ui/c) = 2 
waveguide. The quarter wave resonance condition is how- 
ever only fulfilled for the fundamental mode. The total 

reflection coefficient |r*° tal | will therefore oscillate as a 
function of the number of periods, depending on the in- 
terference conditions for the energy transmitted into the 
high number of higher order modes. 

For grating structures of increasing lateral width, the 
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FIG. 9. The reflection coefficient r*i as a function of the 
number of periods in three GaAs/Alo.3Gao.7As Bragg grat- 
ings. The central waveguide has width a{ui/c) = 1. The 
middle curve displays the reflection properties of a grating of 
lateral width a(u/c) — 2. The reflection coefficient for a sub- 
strate grating and a grating of the same width as the central 
waveguide is also shown for reference (taken from Fig. [5] and 
Fig.0. 



wavevector separation between neighboring modes will 
decrease. The fluctuations will therefore be smoothed 
out in the limit of very wide grating structures. 



VI. RESULTS, NANOWIRE STRUCTURE 

In this section we consider the effect of Bragg reflec- 
tors at the end facet of a semiconductor nanowire. A 
GaAs nanowire will be used as an example, with a Bragg 
grating consisting of GaAs/Al 3 Ga .7As. The nanowires 
have a hexagonal cross-section. The lateral size of the 
waveguide is described using the effective radius p e ff, de- 
fined such that a hexagon with effective radius p e g has 
the same area as a circle with radius p e s- 

Substrate Bragg gratings were shown in the previ- 
ous section to be inefficient for highly diffractive pla- 
nar GaAs/AlGaAs waveguides. This was because only 
a small portion of the fundamental mode is transmit- 
ted into propagating plane waves with high reflectivity 
when diffraction effects arc significant. In a 2D waveg- 
uide there will in general be coupling between modes of 
various polarization, thus both the TE and TM plane 
waves of the substrate have to be taken into considera- 
tion. Fig. [TU] shows the reflection spectrum for TE and 
TM plane waves incident towards a substrate Bragg grat- 
ing. For both polarizations there exists a region with 
reduced reflectivity for transverse wavevectors between 




1 2 

Transverse wavevector. 



FIG. 10. Reflection coefficients, \m\ , for plane wave com- 
ponents in an Alo.3Gao.7As substrate. The plane waves are 
incident toward a 100 period GaAs/Alo.3Gao.7As Bragg grat- 
ing terminated by air. 



0.75(w/c) and (ui/c). This is a strong indication of the 
limited effect of Bragg reflectors in the substrate. A fur- 
ther study of Bragg gratings in the substrate is therefore 
omitted. Note however that due to the Brewster effect 
one might achieve very low reflectivity for TM polarized 
light in the absence of a Bragg grating. A Bragg grating 
may therefore help somewhat in those cases. 

Bragg gratings within the waveguide were seen to be 
promising in the planar case. The analysis is therefore 
extended to find the reflection properties of a hexago- 
nal nanowire using such a grating. The two alternating 
materials are again taken to be Alo.3Gao.7As and GaAs. 
The modes of the hexagonal waveguides were found using 
Comsol Multiphysics™. When extending from a ID to 
a 2D analysis, there is a large increase in computational 
resources. We have therefore chosen to limit the 2D cal- 
culations to a more qualitative analysis, i.e. the simula- 
tion parameters are such that some inaccuracy should be 
expected in the results; we limit the number of modes in- 
cluded in the calculation and reduce the computational 
cell size. Two waveguide widths have been considered; 
p e f[Uj/c = 1 and p e gcj/c = 10. In the simulation we used 
L = 25 and 7 = 0.1; /3u m = ( w /c) for p c sui/c = 10 and 
Aim=0 for peguj/c = 1. Fig. [TTI displays the reflectivity 
of the fundamental mode as a function of the number 
of layers in the grating. The behavior is very similar 
to what was seen for the planar case. The grating does 
yield a significant increase in reflectivity. For the largest 
waveguide, the reflectivity of the fundamental mode as a 
function of the number of periods is in fact almost identi- 
cal to what was seen for the corresponding planar waveg- 
uide. For the smaller waveguide, where diffraction effects 
are more pronounced, the reflectivity converges towards 
a value around 0.93. The deviation from the convergence 
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FIG. 11. The reflection coefficient |ri° tal | 2 as a function of the 
number of periods in the GaAs/Alo.3Gao.7As Bragg grating, 
for three waveguides of various width, a. 



limit of the corresponding planar structure is within the 
uncertainty due to the limited simulation parameters. 

The maximum width for single mode operation in the 
(ZB) GaAs nanowire is around p e ff(w/c)=0.7. This im- 
plies an efficient radius of 97 nm when the excitation 
light is at the lasing wavelength A=870 nm. The cor- 
responding length of one period in the grating would be 
237 nm, and a grating of 20 periods would thus be 4.7 um 
long. Such long nanowire gratings might be challenging 
to achieve with today's technology. For practical pur- 
poses one is thus limited to grow much shorter gratings. 
Higher reflectivity for shorter gratings may be achieved 
by increasing the index contrast, as is e.g. clearly seen in 
Fig. O An increase in the aluminum composition up to 
x = 0.7 would yield a refractive index of 3.15. In a pla- 
nar structure with a(uj/c) = 1, such a grating would be 
capable of achieving a reflectivity of 0.9 after 7 periods, 
compared to 16 periods for x = 0.3. A similar increase is 
to be expected for nanowire structures. It might be pos- 
sible to also perform wet etching of the AlGaAs layers, as 
has been successfully done for VCSEL o 16 i 17 . Wet etch- 
ing of AlGaAs with high aluminum content will increase 
the refractive index contrast further, the refractive index 
of the oxidized AlGaAs layer is around 1.6^. A reflec- 
tor consisting of such oxidized layers would however be 
non-conductive. Wet etching would thus be challenging 
for electrically driven nanowire lasers, as the end facets 
cannot be used for current injection. 



VII. CONCLUSION 

A semi-analytical model has been used to analyze the 
reflection properties of Bragg reflectors to increase the 
end facet reflectivity of diffractive waveguides. Such grat- 
ing are promising to enable high quality nanowire laser 



cavities. We have considered a geometry with the grat- 
ing within the waveguide/nanowire itself and a geometry 
with a substrate grating. The substrate grating has the 
advantage that the composition and thickness are more 
easily controlled, compared to the grating within the 
waveguide. For diffractive waveguides they were however 
found to yield a surprisingly small reflectivity. On the 
other hand, using the geometry with the grating within 
the waveguide, one could obtain near unity reflectivity 
even for extremely small waveguides, where only a small 
fraction of the field is within the waveguide. This would 
however imply using a high number of periods. 

The semi-analytical model enables us to understand 
the mechanisms governing the efficiency of reflection 
gratings in connection to diffractive waveguides. The 
model is however not able to give very exact results for 
two dimensional highly diffractive waveguides, unless a 
high number of radiation modes and evanescent modes 
are included. 

The structure with the grating within the waveguide 
is clearly seen to be the most promising when diffrac- 
tion is significant. For GaAs waveguides terminated 
by such GaAsZAlo.3Gao.7As gratings, one obtains max- 
imum reflectivity after approximately 40 periods of the 
GaAs/Alo.3Gao.7As grating, both for the planar waveg- 
uides and the nanowire waveguides. This is a high num- 
ber that is not easily achieved with today's technology. 
To reduce the number of periods in the grating, one might 
increase the refractive index contrast. For the example 
considered here this could be achieved by increasing the 
aluminum composition in the AlGaAs layers, possibly in 
combination with wet etching. This will create a steeper 
increase in reflection as a function of the number of pe- 
riods, but will also somewhat reduce the maximum ob- 
tainable reflectivity. 
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Appendix A: Reflection and transmission at boundary 
between waveguides 

This section describes the procedure for calculating 
the reflection and transmission matrices at a boundary 
between two waveguides. The situation is sketched in 
Fig. HU Let waveguide b be in the half-space z < 0, and 
waveguide c be in the half-space z > 0. The electro- 
magnetic field is discretized at both sides of the interface 
using the waveguide modes. Note that the modes can be 
divided into two parts; the discrete bound modes, and the 
continuous radiation modes. An artificial boundary con- 
dition e.g. periodic or metallic has to be applied to both 
half-spaces in order to fully discrctizc the modal spec- 



11 



n , b 

Cl I 


, A cl 


tc,b 

n b J 

waveguide b 


>■ n c 

CO 

tb,c 

^ ^ waveguide c 







{e,.hj 



FIG. 12. Geometry for reflection between two waveguides & 
and c. 



trum. The modal field of waveguide b is described by the 
set {ei, hj}, where ei = e%(x, y) and hi = hi (a;, y) are the 
electric and magnetic fields, respectively, of mode i. Sim- 
ilarly, the modal field of waveguide c is described by the 
set {ei,hi}. The reflection and transmission matrices at 
the interface can now be found using the boundary con- 
ditions. Assume that mode i of waveguide b is incident 
from the left, the continuity of the transverse electric field 
yields 



»(*) 



E^H t} 



_ ,cbg(t) 



(Al) 



valid for all x and y. Here r^ b is the reflection coefficient 
from mode i (z < 0) to mode j (z < 0), when the mode 
is incident from waveguide b toward waveguide c, and 
is the corresponding transmission coefficient from mode 
i {z < 0) to mode g (z > 0). The superscript (t) stands 
for the transverse component (x and y components) of 
the vector. The boundary condition for the transversal 
magnetic field is similarly 



,(*) 



E 



r cfc h (t) = V < cb h^ 



(A2) 



Eqs. (|A1|) and (|A2[) can be combined as follows. Take the 

vector product between (|A1[) and Hm* {x,y), and inte- 
grate over the unit cell. Similarly, take the vector product 
between Sm*(x, y) and (|A2[) , and integrate over the unit 
cell. Here, { / H m ,£m} belong to mode m of an ambient 
dummy medium with refractive index n a , according to 
Eq. ©. This yields 



,cb 



1 



mq ■ 



(A3a) 
(A3b) 



Here, ^ m i, *&mi and f m i, <& m i are inner products, as 
given by Eq. (fT2|) . between the free space modes (O and 
the modes of the waveguides, i.e. {ei,hi} and {ei,hi}, 
respectively. Let t ab = [t mi ] ab and r ab = [r ji ] ah be the 
shorthand notation for the reflection and transmission 
matrices for light incident from waveguide b (with modal 
fields {ei,hi}) towards the homogeneous ambient, and 
similarly for waveguide c (with modal fields {ei,hi}). 
From (fTTI) we have: 

2t ab = ([*rm] + I*™]) - - [*mi})r ab (A4a) 

2t ac = ([* mi ] + [$ mi ]) - - r ac (A4b) 

(W + r Qb = ([$ roi ] - [*rm]) (A4c) 

([**.<] + r oc = ~ [*mi]) ■ (A4d) 

We eliminate ([* mi ] - and ([# mi ] - [$„„]); 



2t° 6 (/-r a( Y o6 ) + = ([$ mi ] + [*„,,]) 



(A5a) 

2t ac (/ _ r ac r ac )+ = (^.J + _ (A5b) 

Eqs. (|A"3a| and (|A"3b")) can now be rewritten 

([**.<] + [*mi]) {I - r ab r cb ) = (l$mi) + $mi]) t ch 

(A6a) 

([*mi] + [*mi]) - 7- C ") = ([*mi] + [*mi]) r 4 "^. 

(A6b) 



We solve (|A"6a|) and (jA"6b)) for r ch and i cb ; 



c& s~icb 



G c0 J-r a V 



(A7a) 



c6 { ficb ^ac/~icb^ab\ Jr fQcb^ab ^ac£icb\ (A71o) 



(G c0 - r ac G c V 
Here we introduced 



G cb = ([$ mi ] + [<F mi ]) + ([$ mi ] + [*„,,]) (A8a) 



= (/ - r ac r ac ) (t ac ) + t ab (I - r ab r ab ) 



(A8b) 



The opposite reflection and transmission, r bc and t bc , 
can be found in the exact same way. The result is im- 
mediately available by interchanging the indices b and c, 
i.e. 



be /~ibc 



q c j _ r ac r , 



ac^bc\ 



(A9a) 



r bc = ( G bc _ r ab G bc r ac^+ (g&c^c _ ^fcg&c) ^ ( Agb ) 

where 

G bc = (G cb ) + . (AlOa) 



Using Eqs. (|A9|) and Eqs. (|AT[) . in combination with the 
recursive relations ([5]), the scattering properties of e.g. a 
complete Bragg structure can be calculated. 
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